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Data on \ow-px hadronic spectra are widely regarded as evidence of a hydrody- 
namic expansion in nucleus-nucleus collisions. In this interpretation, different hadron 
species emerge from a common medium that has built up a strong collective velocity 
field. Here, we show that the existence of a collective flow field implies characteristic 
modifications of high-py parton fragmentation. We generalize the formalism of par- 
ton energy loss to the case of flow- induced, oriented momentum transfer. We also 
discuss how to embed this calculation in hydrodynamic simulations. Flow effects are 
found to result generically in characteristic asymmetries in the r/ x ^-plane of jet 
energy distributions and of multiplicity distributions associated to high-p-p trigger 
particles. But collective flow also contributes to the medium-induced suppression of 
single inclusive high-p^ hadron spectra. In particular, we find that low-pr elliptic 
flow can induce a sizeable additional contribution to the high-py azimuthal asymme- 
try by selective elimination of those hard partons which propagate with significant 
inclination against the flow field. This reduces at least partially the recently observed 
problem that models of parton energy loss tend to underpredict the large azimuthal 
asymmetry V2 of high-pr hadronic spectra in semi-peripheral Au+Au collisions. 



1. Introduction 



What happens if a hard process, such as the production of high-.&r jets, is 
embedded in a dense nuclear environment created e.g. in a nucleus-nucleus 
collision at RHIC or at the LHC ? While parton-parton interactions at high 
virtuality Q 2 ^> A.q CD occur on too short time and length scales to be af- 
fected by the typical modes in the medium, the parton showers associated to 
the incoming and outgoing state can interact with the medium [1-6]. This is 
expected to result in an energy degradation of the leading parton [1-6], in a 
transverse momentum broadening of the parton shower [7-10], and in an en- 
hanced and softened multiplicity distribution of the hadronic final state [10]. 
Measurements in Au+Au collisions at RHIC support this picture by the ob- 
served suppression of leading hadron spectra [11-16] and leading back-to-back 
correlations [17], as well as the medium-modified "jet-like" properties of parti- 
cle production associated to high-p T trigger particles [18-20]. The analysis of 
these "jet quenching" observables has become one of the most active and most 
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diverse research fields in ultrarelativistic nucleus-nucleus collisions, mainly be- 
cause the pattern of medium-induced partonic energy loss is expected to allow 
for a detailed characterization of the properties of the produced dense medium. 

Parton energy loss is known to be sensitive to the total in-medium path- 
length and to the average squared transverse momentum transferred from the 
medium to the hard parton [21-25]. In recent phenomenological studies, the 
latter quantity has been parameterized by the BDMPS transport coefficient 
q [10,26] or by physically equivalent model-dependent quantities such as twist- 
4 multiple scattering matrix elements [27], the medium opacity or the number 
of initially produced gluons per unit rapidity [28]. These model parameters 
can be related to the energy density of the produced matter [10,29,30]: 

?(0 = cc s/4 (0- (1-1) 



Here, c is a medium-dependent constant (c ~ 2 for the case of an ideal quark 
gluon plasma [26,29]), and (1.1) also holds for a time-dependent energy density 
as indicated by the explicit ^-dependence. Recently, several phenomenological 
analysis have used the medium modifications of high-p^ hadron production to 
extract information about the energy density attained in nucleus-nucleus col- 
lisions at RHIC [24,26,27,31]. These models also account successfully for the 
centrality dependence of the suppression pattern [27,31,32], and the reduction 
of leading back-to-back correlation [27,31,32]. However, they tend to under- 
predict [31,32] the elliptic flow V2{pt) at high transverse momentum which is 
thought to originate from parton energy loss in an azimuthally asymmetric 
geometry [28,33]. 

In nucleus-nucleus collisions at RHIC and at the CERN SPS, there is strong 
experimental evidence that the produced medium is - if at all - only locally 
equilibrated and is thus characterized only locally by its energy density. Mea- 
surements of 1ow-pt inclusive hadron spectra [34,35] and their azimuthal asym- 
metry [36-38] support the picture that different hadron species emerge from a 
common medium which has built up a strong collective velocity field [39-43]. 
These measurements are broadly consistent with calculations based on ideal 
hydrodynamics [41-45], in which the dynamic behavior of the produced QCD 
matter is fully specified by its equation of state p = p(e,T,/i B ) which enters 
the energy momentum tensor 

T^(x) = (e + p) u^u v -p . (1.2) 



For the case of a longitudinal Bjorken-type flow field = /Vl — P' 2 -, 

— * 

(3 — (3 z, the longitudinal component of the energy momentum tensor increases 
from T zz = p to T zz = p + Ap, where Ap — (e + p)u z u z = 4p/3 2 /(l — (3 2 ) 
for the equation of state of an ideal gas, e = 3 p. For a rapidity difference 



2 



r] = 0.5, 1.0, 1.5 between the rest frame which is longitudinally comoving with 
the jet, and the rest frame of the medium, the component T zz of the energy- 
momentum tensor "seen" by the hard parton is thus increased by a factor 
1, 5, 18, respectively. It is thus reasonable to assume that the momentum trans- 
fer from the medium to a test particle such as a hard parton does not depend 
solely on the local energy density e, but rather on the energy momentum tensor 
(1.2) which involves a significant directed collective flow field u^{x) [46]. 



Vacuum Static medium: Flowing medium: 

(reference) Broadening Anisotropic shape 




Fig. 1. Sketch of the expected energy or multiplicity distribution of a jet frag- 
menting i) in the vacuum (left), ii) in a medium which is longitudinally comoving 
with the rest frame of the jet (middle) and iii) in a medium which is longitudinally 
boosted with respect to the rest frame of the jet (right). 

Figure 1 sketches the qualitative picture first advocated in Ref. [46]: A jet 
which fragments inside a medium is known to broaden its shape and to soften 
its multiplicity distribution. But if the medium exhibits a collective motion, 
then a smaller local energy density is sufficient for the same net momentum 
transfer to the hard parton and thus for the same medium-induced parton 
energy loss. Moreover, the directed momentum transfer can be expected to 
break the rotational symmetry of the jet shape in the rj x </>-plane. In this work, 
we give a detailed description of the formalism incorporating these effects and 
we explore observable consequences of the resulting interplay of oriented and 
random momentum transfer to a hard parton. 

For each jet, rotational symmetry in the rj x 0-plane is broken even in the ab- 
sence of a medium, mainly for two reasons: Statistically, any finite multiplicity 
distribution of a rotationally symmetric sample breaks the rotational symme- 
try. If this were the only source of symmetry breaking, one could search for 
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medium- induced asymmetries in realigned jet samples, similar to the analysis 
of elliptic flow in realigned event samples [47, 48] . In addition, however, the 
fcr-ordering of the DGLAP parton shower implies that the first parton split- 
ting in the shower contains significantly more transverse momentum than the 
second, thus leading to a dynamical asymmetry in the rj x 0-plane. Both effects 
lead to a symmetry breaking in a random direction in the 77 x 0-plane - thus 
rotational symmetry is restored in sufficiently large jet samples. To search for 
symmetry breaking effects caused by collective motion in 77 x 0-distributions of 
jet energy and jet multiplicity, it is thus important to control experimentally 
the direction of this collective motion. Based on these arguments, we foresee 
two classes of applications for our calculations: 




Fig. 2. Schematic view of two scenarios in which jets interact with collective 
flow fields: a) If the hard parton is not produced in the Lorentz frame longitudinally 
comoving with the medium, or if the longitudinal collective flow does not show 
Bjorken scaling, then the parton interacts with a flow component parallel to the 
beam, b) On its propagation in the transverse direction, hard partons generically 
test transverse flow components, except for the special trajectories which are parallel 
to the flow field. 

First, in general, a hard parton needs not be produced in the Lorentz frame 
which is longitudinally comoving with the medium; and even if it is produced 
in the longitudinally comoving frame, it will in general not stay in this frame 
during the entire time evolution of the medium. This is so since the hard 
parton moves - like any effectively massless particle - on a straight light-like 
line in the (z, t)-diagram, whereas the collective flow field is expected to show 
significant deviations [49, 50] from Bjorken expansion and will thus intersect 
this straight line. In such cases, the collective component of the momentum 
transfer to the hard parton is directed along the beam axis. Hence, averaged 
samples of medium-modified jet shapes and jet multiplicities can be expected 
to show an asymmetry which is preferentially oriented along the beam direc- 
tion in the 77 x 0-plane. [At mid-rapidity, the jet sample must be symmetric 
with respect to the 77 — > —77 mirror symmetry, but - in general - it will not 
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be rotationally symmetric in the rj x 0-plane.] In the case of a significant 
transverse collective flow, the analogous argument implies the occurrence of 
jet asymmetries preferentially oriented along the ^-direction. In general, the 
size and orientation of the jet asymmetry depends on how hard processes are 
spatially distributed in the dynamical expansion scenario. 

Second, flow effects manifest themselves not only in the azimuthal asymme- 
tries of jet observables, but also in inclusive high-p^ hadron spectra. In the 
presence of a flow field, a smaller local energy density is sufficient for the same 
net momentum transfer to the hard parton and thus for the same medium- 
induced parton energy loss. This is relevant for the interpretation of the nu- 
clear modification factor in terms of a local energy density. Moreover, hard 
partons propagating parallel to a flow field can be expected to suffer less mo- 
mentum transfer and hence less parton energy loss than those traveling along 
non-parallel trajectories. In this way, low-p^ collective flow can induce high- 
p T azimuthal asymmetry by selective elimination of those hard partons which 
propagate with significant inclination against the flow field. We show that this 
can yield to a sizeable additional contribution to high-p T t> 2 . 

This paper is organized as follows: In section 2, we introduce the formal- 
ism in which the effects of anisotropic momentum transfer on parton energy 
loss are calculated. In section 3, we calculate the induced asymmetry of the 
medium- dependent gluon energy distribution and in section 4 we analyze the 
resulting anisotropic jet energy distribution. Up to this point, the medium 
will be characterized schematically by the momentum scale /x which deter- 
mines the random momentum transfer to the hard parton, and the vector q 
which specifies the oriented momentum transfer. The ratio qo/fi indicates the 
relative strength of collective flow. In section 5, we then discuss how to embed 
this calculation in a dynamical expansion scenario and we estimate the flow- 
induced parton energy loss contribution to high-pr V2- Finally, we summarize 
the main conclusions. 



2. The formalism 



The starting point of our calculation is the energy distribution of gluons into 
which the initially produced parent parton fragments, 

dI tot E T -AE T dl™ c dl mcd . . 

uj = uj hw . (2.1) 

dudk E T dudk dudk 

Here, E T denotes the total energy of the hard parton and AE T = 
J duj dkuj jf^K- is that part of the total energy which is redistributed by 
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medium- induced radiation. The factor (E T — AE T ) / ' E T ensures energy- 
momentum conservation. 

In the absence of a nuclear environment, the parton fragments according to 
the distribution I tot = I vac . In the medium, the parent parton radiates ad- 
ditional gluons due to medium-induced multiple scattering. This medium- 
induced gluon radiation has been calculated to leading order in l/E, resum- 
ming an arbitrary number of scattering centers [2-5]. It depends to leading 
order in l/E on the in-medium pathlength L and on the average squared 
transverse momentum transferred to the hard parton per unit pathlength. 
The latter property of the medium is parameterized differently in different ap- 
proaches, e.g. in terms of the BDMPS transport coefficient q [2-4], or as the 
product of the longitudinal density of scattering centers Hq along the parton 
trajectory, and their typical momentum transfer fj 2 [4,5]. These parameteri- 
zations are known to lead to equivalent results for the medium-induced gluon 
radiation [9]. In what follows, we use the single-hard scattering approximation, 
in which the effects of multiple scattering are characterized by the number of 
effective scattering centers Uq L times the radiation off a single scattering cen- 
ter. The elastic scattering cross section is modeled in terms of Debye-screened 
Yukawa-potentials 

The medium-induced gluon radiation is given by [4,5,9] 

dl™« _ a s 4C R n r k . q -L^ + sin(L^) 

U d^-j2^—^- jW^I I? [(k + q ) V 2c] 2 • 

(2.3) 



This radiation spectrum is for a time-independent homogeneous medium. 
However, it also applies to a time- dependent density of scattering centers if 
the density n is replaced by an appropriate time-average (see Ref. [9,51] and 
section 5 A below). 

If a flow component is directed orthogonally to the parton trajectory, then 
the momentum transfer from the medium is anisotropic. We denote by 
q = (q , q l Q ) the directed momentum transfer to the hard parton which is 
parallel to the spatial components of the flow field u^{x). Here and in the 
following, transverse vectors lie in the plane orthogonal to the trajectory of 
the hard parton while longitudinal components are parallel to this trajectory. 
In the high-energy limit, momentum transfers parallel to the hard parton are 
negligible. Thus, the effect of collective flow on the medium-induced radiation 
(2.3) can be accounted for by using an anisotropic scattering potential 



6 



\a{q) 



7T 



(q-qo) +^ 



2 ' 



(2.4) 



The parameters /x and |qo| characterize the strength of the random and di- 
rected momentum transfers from the medium to the hard test particle, re- 
spectively. The component q l is parallel to the parton trajectory and does not 
enter our calculation. We work in radial coordinates, 

dq = q dq dip , d\t = kdkda, (2-5) 



where a denotes the angle between the transverse momenta k and qo- The <p- 
integration in (2.3) can then be done analytically. We always work in the frame 
longitudinally comoving with the hard parton in which the parton propagates 
orthogonal to the beam direction. 
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Fig. 3. Definition of kinematic variables of a gluon emitted inside a jet cone. 
Variables are defined in the Lorentz frame which is longitudinally comoving with 
the jet. 



We now consider more explicitly the case of a longitudinal flow component 
parallel to the beam direction, see Fig. 3. To express the radiation spectrum 
as a function of pseudorapidity rj and azimuthal angle <fi with respect to the 
center of the jet at rj = = 0, we write the pseudorapidity of an emitted gluon 
as rj = — In tan | where (3 is the angle between the gluon momentum and the 
beam axis. We have 



tan (3 



UJ 



k cos a \ 



cos 2 a . 



tan< 



k sin a 
v^ 2 - k 2 



(2.6) 



7 



Inversion leads to 



cos a 



sinh?7 



k \J cosh 2 rj — cos 2 <p 



\J cosh 2 1] — cos 2 <p 



to 



cosh?? 



(2.7) 



The Jacobian for the transformation to jet observables 77, reads 

COS(' 



kdkda = uo 



cosh 77 



G?77 d(p . 



(2.8) 



Our final expression is 



u dl™ d cos0 a^C R 4n /x 2 /-qo ^a-- s jn^- 

dudrjdcj) cosh 3 ?? 7r 2 k 2 io 1 q 4 

k 2 (g 2 + /i 2 + (q + k) 2 ) - 2q 2 k ■ (q + k) 



(q 2 +/i 2 -(q + k) 2 ) +4/i 2 (q + k) 5 



3/2 ' 



(2.9) 



where k and the angle a are given by (2.7) in terms of 77 and <p. The case of a 
transverse flow component is obtained by obvious rotations and redefinitions 
of the vectors in Fig. 3. 



3. Transverse momentum dependence of the medium-induced gluon 
radiation in the presence of flow 



In this section, we discuss the generic properties of the medium-induced gluon 
radiation (2.3) in the presence of collective flow. To this end, it is convenient 
to change to dimensionless variables 

K=\k\/n, q=\q\/n, 7 = u c /uj, tu c = ^/j 2 L . (3.1) 



We transform (2.3) to radial coordinates dq = fi 2 qdqdip and dk = fj, 2 R,dR,da, 
where a denotes the angle between the transverse momenta k and q . Doing 
the (^-integration in (2.3), we find 



LO 



dl™ d _a s C Rn r 2 f ~ \ sin 7g 2 



2n L J dq 



doj k dn da n 2 J k 2 q 4 

R 2 (q 2 + 1 + (£0 + r2 + 2 <?o« cos a)) - 2q 2 (R 2 + Rq cos a) 

*7 3/2 ■ \^ ) 

(q 2 + 1 — (<7q + k 2 + 2q Rcosa)) + 4(q 2 + k 2 + 2q Rcosa) 
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Fig. 4. Plot of the medium-induced gluon energy distribution (3.2) as a function 
of the angle a and rescaled transverse momentum R 2 = k 2 /^ 2 for different values 
of the rescaled gluon energy ^/x 2 L/u and the collective flow strength q$. Plots are 
for jj, = 2 GeV and L = 6 fm. 



In Fig. 4, we plot the medium-induced energy distribution as a function of 
R 2 and a for fixed ratios of 7 = lj c /uj. Here, a = denotes the direction 
of the collective flow vector qo. The figure 4 shows clearly that more energy 
is deposited in the direction of the collective flow vector qo. For the same 
reason, the energy distribution is depleted in the direction opposite to qo, i.e. 
for a ~ 7T. For the medium-induced component plotted in Fig. 4, this shows 
up as a negative contribution, while the total energy distribution (2.1) stays, 
of course, positive. 

For non-zero values of the flow vector q , the triple differential gluon distri- 
bution (3.2) has a singular behavior for R — > 0, 

dl™ d / +00 for -!<«<!, 
s^o duondnda \ -00 tor 2<a<y. 

Figure 4 displays only finite values of R, but the limit (3.3) is consistent 
with the small-/? behavior seen in Fig. 4. This singularity is unphysical. It 
stems from the fact that the formalism leading to (2.3) calculates medium- 
modifications to a perturbative parton splitting oc ^ without regularization 
of this collinear divergence. At R — 0, the anisotropic flow field shifts part of 
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this singularity as a positive contribution to the half plane — ~ < a < ~ while 
depleting the region | < a < In the x 0-plane, (2.9) integrated over 
energy shows the same singular structure, namely a positive divergence for 
?7 — ^ + and a negative divergence for rj — > 0~. These two divergences cancel 
each other if integrated over an arbitrary small neighborhood around R = 
(i.e. around rj = cf) = 0). They represent a very small contribution to the total 
jet energy (see the discussion below). 




Fig. 5. The gluon energy distribution (3.2) for the same parameters as in Fig. 4, 
but averaged over the cases qo and — qo- 

In Fig. 5, we plot the gluon energy distribution (3.2) for a sample of medium- 
modified jet shapes for which the orientation ±qo of the collective flow is 
known but the directions qo and — qo are equally likely. The figure clearly 
indicates that the singular behavior for R — > disappears. For fixed values 
of 7 = — , the energy distribution has a finite constant value for R — > 0. By 
construction, the obtained distribution is symmetric around a — |, but it 
shows a marked angular dependence. 



4. Collective flow leads to anisotropic jet energy distributions 

In this section, we discuss the medium-modification of jet energy distributions. 
Our starting point is the gluon energy distribution (2.1). We constrain the 
vacuum contribution J vac of this spectrum by data on the energy fraction of 



10 



a jet contained in a subcone of radius R = y/r/ 2 + (ft 2 , 
„ f&\ - E T (R) 



Njcts E T (R = 1) 

= du; c?k ( — - — i? ] -j- — — . (4.1) 

E T J J \u J dudk K J 



For the jet shape p vac (R), we use the parameterization [52] of the Fermilab 
DO Collaboration determined for jets in the range ~ 50 < E t < 170 GeV and 
opening cones 0.1 < R < 1.0, 



p V£) (R) =AR 31 + B R 06 + C R"- b + D R 7 + E i? uy , (4.2) 
where 

A(E T ) = -3.47 + 0.85 • 10~ 2 £ T - 0.25 • 10" 4 £| , 
D(E T ) = 3.30 - 0.77 • 10~ 2 £ T + 0.22 • 10" 4 £| , 

5 = 9.75, C = -8.32, £ = -0.30. (4.3) 

At very small values of R, this parameterization turns negative. We cure this 
unphysical behavior by matching (4.2) with a third order polynomial which 
smoothly interpolates to p vac (R = 0) = 0, 

„ (R )-I P^cKR) for R>0M, 

P^ R )-\ aR 2 + bR* for i?<0.04. (44) 



The parameters a and b in (4.4) are fixed by requiring that p vac (R) and its 
first derivative are continuous. 

The medium-induced part of the jet energy distribution is determined from 
(2.9), 

dE med fE dI mcd 

d UUJ TTIZTJI ■ ( 4 - 5 ) 



dr) d<p Jo duo dr] 

It contains a fraction of the available jet energy, 

For the vacuum contribution, the corresponding distribution is defined by the 
vacuum jet shape (4.4), 
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m = {Et ~ AEt) ^ik ■ < 47 > 

Here, the prefactor (E T — AE T ) ensures energy conservation. 
A. Harmonic expansion of jet energy distribution 



In this section, we characterize the flow-induced asymmetry of jet energy dis- 
tributions in terms of a harmonic analysis in the rj x 0-plane. We introduce 
radial coordinates in this plane, 

R= vV + 2 , (4.8) 
a' = arctan ( . (4.9) 



The flow field points in the direction a' = 0. We calculate now the jet energy 
distribution in the 77 x 0-plane, 

dE T dE T dE™ c dE™ d 

T - T 1 T (4.10) 



RdRda' drjdcf) drjdcf) drjd(f) 

In the absence of flow, the radiation spectrum (3.2) is rotationally symmetric 
in the coordinates R (or k) and a. However, the energy distribution (4.10) 
is elongated in the ^-direction due to the Jacobian (2.8) in the coordinate 
transform. In general, this reduces the effect of 77-broadening due to longitu- 
dinal flow, and it can be corrected for analytically. However, this asymmetry 
is rather small (< 10%) for small cone sizes (R < 0.3), and can be neglected 
safely in the following discussion. To analyze the asymmetries of the jet energy 
distribution (4.10) in the 77 x 0-plane, we use a harmonic expansion 

dE 00 

E {0) (R) + 2j2E in) (R) cos (na') . (4-11) 



RdRda' 



The coefficients proportional to sin (na') cancel since the jet energy distribu- 
tion (4.10) is by construction symmetric with respect to a' — > —a'. 
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Fig. 6. The harmonic coefficients which characterize the asymmetry of the jet 
energy distribution of a 100 GeV jet in the 77 x 0-plane. Parameters given in the 
text. 



For illustration, we now study the case of a jet of total energy Et = 100 GeV 
traversing a medium in which AEt = 23 GeV were redistributed in phase 
space due to medium effects. This parton energy loss is obtained e.g. for an 
in-medium pathlength of nuclear size (L — 6 fm), a momentum transfer per 
scattering center of /i — 1 GeV and a collective flow effect of the same size 
go = /i, with an effective coupling constant in (2.3) fixed to rioLa s CR = 1. 
This corresponds to a gluon jet (Cr = Ca = 3) with a reasonable perturbative 
coupling a s = 1/3 and an opacity no L = 1. Changes in the coupling can be 
absorbed in a redefinition of the density of scattering centers. Recent model 
comparisons with RHIC data support the picture of an opaque medium which 
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may result in significantly larger medium- induced modifications than the ones 
modeled here. However, to illustrate the sensitivity of jet shape measurements, 
we prefer to work with a sizeable but relatively small effect. 

In Fig. 6, we show the first 6 harmonic coefficients extracted for the jet energy 
distribution (4.11). Since the vacuum term (4.4) is rotationally symmetric, it 
contributes only to E^(R). The total jet energy is obtained by integrating 
this zeroth moment over R. The shape of the medium-induced part of E^ (R) 
is negative at small R which indicates the medium-induced depletion of the jet 
energy in this region of phase space. This is the result of multiple scattering 
which broadens the jet energy distribution by moving a fraction of the total jet 
energy to larger values of R. The higher moments E^ n \R), n > 1, contain in- 
formation about asymmetries in the energy distribution but do not contribute 
to the total jet energy. For the case of averaged jet samples over opposite 
flow directions qo and — qo, the odd harmonic moments vanish while the even 
ones stay the same. The absolute size of the harmonic moments decreases by 
approximately one order of magnitude per moment n. This indicates that the 
first and second moments are sufficient to characterize the asymmetries of the 
jet shape. We note that to reconstruct from the moments E^(R), E^(R), 
E^ 2 \R) the jet shape, one has to integrate over R dR - this tames significantly 
the large values of these moments in the region R ~ 0. We expect that most 
of the experimentally accessible structures lies in the range 0.05 < R < 0.3. 

The technical advantage of the harmonic expansion (4.11) is that the unphysi- 
cal singularities of the vacuum contribution (4.7) and the medium contribution 
(4.5) at R = can be removed easily. In the harmonic coefficients E^ n \R) of 
(4.11), these singularities appear in the odd moments at R — > 0. The smallest 
value calculated for Fig. 6 is for R = 0.01. For much smaller values of R, we 
find numerically E^\R = 10~ 4 ) ~ 2 • 10 5 and E {1 \R = 10~ 6 ) ~ 2 • 10 7 . For 
the figures presented in this work, we cut off this artificial small- R structure in 
the harmonics and plot the jet energy distribution according to (4.11). As ex- 
plained above, the presence of these singularities indicates that our formalism 
becomes unreliable in the collinear region of very small R. For the arguments 
in this paper, this is not a problem since only a small amount (< 5% for 
R < 0.01) of the total jet energy lies inside this small phase space region. 

B. Profile and displacement of jet energy distribution 

In general, collective flow shifts the calorimetric jet center and distorts the 
shape of the jet energy distribution. Any calorimetric jet finding algorithm 
can be expected to center the cone around the medium-displaced jet center. 
Fig. 7 shows this displacement Arj as a function of the collective flow strength, 
calculated by averaging the jet energy distribution (4.10) over a central cone of 
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size R < 0.3. For a fixed average momentum transfer fj, per scattering center, 
the displacement grows approximately linearly with the directed momentum 
transfer go- Also, At] grows approximately linear with the average momentum 
transfer \x for a fixed ratio qo/fj,. The overall size of the displacement is rather 
small: a displacement of size Arj = 0.1 results only for rather large parameter 
values (e.g. q% = 4 /i 2 = 8 GeV 2 ) which correspond to a large medium-induced 
average energy loss (AE T ~ 62 GeV). This is consistent with the picture 
that the most energetic jet fragments are radiated collinear. These energetic 
components which dominate the energy distribution, are shifted very little in 
the 7] x 0-plane even if they pick up a significant transverse momentum. 




2.5 



3.5 
jj 2 (GeV 2 ) 



Fig. 7. The displacement At] of the calorimetric center of the jet cone as a 
function of collective flow strength qo and average momentum transfer fi from the 
medium. 



In Fig. 8, we show different one- dimensional cuts through the jet energy distri- 
bution (4.10). These cuts go through the displaced jet center and are quantified 
by the radial coordinates Rd and of the displaced center. Along the beam 
direction (a^ = 0), the jet energy distribution is shifted with the flow field, see 
Fig. 8. The medium-induced part of the jet energy distribution takes negative 
values in the region of phase space which is depleted due to medium effects. 
We observe in particular a pronounced long tail of the distribution in the di- 
rection of the flow. We attribute this tail to the soft jet fragments which can be 
displaced significantly in rj by a typical momentum transfer from the medium. 
In contrast, in the direction orthogonal to the collective flow, one observes a 
numerically small medium-induced broadening of the jet energy distribution, 
which is not accompanied by a displacement. 
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Fig. 8. The jet energy distribution (4.10), plotted along different cuts in the 
T] x cf>- plane, as indicated. The variables Rd and a.d denote the distance and orien- 
tation with respect to the displaced calorimetric center of the jet cone. Parameters 
are the same as those of Fig. 6. 



For measurements at mid-rapidity, one can access experimentally only the ori- 
entation but not the direction of the collective flow component. Then, each jet 
of the sample is positioned around its calorimetric center which is shifted with 
equal probability in the positive or negative beam direction. We have given nu- 
merical results for this case already in a previous exploratory study [46]. The 
main conclusion is that the jet energy distribution can broaden significantly 
along the orientation of the flow. The analogous conclusion holds for the case 
of a strong transverse collective flow field, where only — > — ^-symmetrized 
samples are measurable. 
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C. Asymmetries in the jet multiplicity distributions 



If the jet energy distribution is sensitive to collective flow, then high-pr par- 
ticle correlations and jet multiplicities should be sensitive too. Based on the 
medium-induced gluon energy distribution (2.7), we can explore this sensitiv- 
ity qualitatively by calculating the number of medium-induced gluons emitted 
with a gluon energy uo > uj cnt . In the 77 x 0-plane, this medium- induced mul- 
tiplicity distribution reads 



To relate equation (4.12) to an experimentally accessible quantity requires a 
hadronization model and faces uncertainties which have been mentioned be- 
fore [9]. In Fig. 9, we plot the medium-induced modification (4.12) of the jet 
multiplicity for the case of a collective flow in the positive beam direction. 
The qualitatively expected effects are illustrated clearly. In the beam direc- 
tion, there is a marked reshuffling from negative to positive rapidities due 
to collective flow effects. In the direction orthogonal to the beam, there is a 
somewhat smaller reshuffling from smaller to larger cone sizes: this is the result 
of multiple scattering in an isotropic medium which leads to a characteristic 
broadening of the jet multiplicity distributions. 

As a result of the eikonal approximation which underlies the calculation of 
the medium-induced gluon energy distribution (2.7), the leading hard par- 
ton does not change its direction due to medium effects. This corresponds to 
the assumption that the leading hadron of the jet is located at r? ~ ~ 0. 
Then, Fig. 9 provides an estimate of the rotational asymmetry of hadron 
production associated to a high-p^ trigger particle. Recently, two-particle cor- 
relations and their possible medium-modifications have been discussed in sev- 
eral approaches [53-57]. Based on the present study, we expect flow-induced 
asymmetries to affect such two-particle correlation measurements in nucleus- 
nucleus collisions. In particular, asymmetries are expected in multiplicity dis- 
tributions associated to high-p^ trigger particles and in leading two-particle 
correlations [46]. 



N((j>, ?7, Wcut) 



dN(u mt ) 
drj d(p 
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Fig. 9. The medium-induced contribution (4.12) to the total jet multiplicity- 
distribution for different cuts in n x </>-plane. Parameters are the same as those of 
Fig.6. 



5. Parton energy loss in dynamical simulations 

A. A proposal to determine parton energy loss from the energy-momentum 
tensor 

In a realistic dynamical scenario of a nucleus-nucleus collision, the produced 
hard parton propagates through a medium of varying spatial and temporal 
energy density e(r, z, £) and varying collective flow u^(r, z,£). Thus, its parton 
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energy loss will depend on the spatial position r of its production point (the 
production time is £ ~ for hard processes), and the orientation n of its 
trajectory, 

m=r + £n. (5.1) 



In the absence of collective flow, numerical studies of the medium-induced 
gluon energy distribution have shown that the medium- induced gluon ra- 
diation for a medium of time- dependent density is equivalent to that of a 
static medium whose density has been rescaled appropriately. This rescaling 
requires the determination of the linearly line-averaged characteristic gluon 
energy along the trajectory r(£) [28,51,58], 

roo 

Jo 



Here, we have expressed the BDMPS transport coefficient (1.1) in terms of the 
local energy density, q(f, £) = ce 3//4 (r, £). In the same way, one can determine 
the time-averaged total transverse momentum squared 

{qL)m\ = / ^ce 3 / 4 (r(0,0, (5.3) 
Jo 



and construct the quotient to c L = 2|g. The probability distribution that an 
additional fraction AE of the parton energy is lost due to medium-induced 
scattering depends on uj c and uj c L. A numerical routine for its calculation is 
publicly available [9]. The characteristic gluon energy (5.2) and momentum 
broadening (5.3) can be related to the model parameters entering the medium- 
induced gluon energy distribution (2.3) via 

u c = X -qL 2 = ^(n L) /i 2 L , (5.4) 
qL = (n L) /i 2 . (5.5) 



Remarkably, the time-averages (5.2) and (5.3) do not require an a priori knowl- 
edge of the in-medium pathlength L. In the case of a time-independent energy 
density of a medium of finite size, e(r) oc 0(|r| — L), one recovers the expres- 
sions for the static case. For further details of how to relate parton energy loss 
in a time- dependent medium to an equivalent time-independent calculation, 
we refer to Ref. [9,51,59]. 

Collective flow is an additional source of momentum transfer to the hard 
parton, and will result in additional parton energy loss. To account for this 
effect, we suggest to replace the energy density in (5.2) by the relevant boosted 
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component of the energy-momentum tensor (1.2). To be more specific, we 
consider the component T n±n± where n± is orthogonal to the trajectory (5.1) 
of the hard parton, 

r^=p(c) + [e + p(e)]^-. (5.6) 



Here, (3± is the spatial component of the collective flow field which is orthog- 
onal to the parton trajectory. In general, all quantities entering (5.6) will de- 
pend on space and time. In the absence of flow effects, f3± — 0, the component 
T n±n± = p determines the pressure and hence it determines via the equation 
of state the energy density e{p) entering (5.2) and (5.3). For finite flow 
our proposal is to use e(T n±n± ) instead of e(p) in evaluating the characteristic 
gluon energy and momentum broadening, 

g = ce 3 / 4 (p) _ > q = ce 3/4 (T n ^). (5.7) 



This is consistent with what is known from analytical estimates and numeri- 
cal studies about the dependence of parton energy loss on momentum transfer 
from the medium. For the determination of jet asymmetries in a dynami- 
cal scenario, relation (5.6), one has to determine the relative strength of the 
random and directed momentum transfers in (2.3). For a viable model, qo/fi 

should increase monotonically with t+ p P ^ y^ff 



2 • 



B. Low-pr elliptic flow induces high-px azimuthal asymmetry 



In general, a hard parton will suffer less energy loss if it propagates on a trajec- 
tory parallel to the flow field. Thus, for the same medium-induced suppression, 
the azimuthal asymmetry at high transverse momentum becomes larger when 
the contribution of the collective flow field is increased. To estimate the size 
of this effect, we consider a simple two-dimensional model. The hard parton 
is produced at an arbitrary position (x ,yo) m the transverse plane according 
to the nuclear overlap. It propagates in its longitudinally comoving rest frame 
in the transverse direction n = (cos <p, sin <p), along the trajectory 

r (0 = (x + £ cos tp, yo + £ sin tp) . (5.8) 



For simplicity, we assume that the longitudinally comoving rest frame of this 
hard parton is the longitudinal rest frame of the medium. Then, there is only a 
transverse but not a longitudinal flow component. For the BDMPS transport 
coefficient which includes collective flow effects, we make the ansatz 
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9(0 = Qnf + ?/|«r(r (0) ' n r 



(5.9) 



Here g/ and q n f stand for the flow and non-flow components to q, the two- 
dimensional vector n T is orthogonal to the trajectory (5.8) and projects out 
the corresponding transverse component of the collective flow field ut(to(£)). 
We discuss now the motivation for this ansatz. In the absence of collective 
flow, q n f defines the time-averaged BDMPS transport coefficient of the dy- 
namically equivalent static scenario, as specified in the discussion of (5.2) 
and (5.3). Thus, the ansatz (5.9) can account for one of the main effects 
of longitudinal expansion, namely the time-dependent decrease of scatter- 
ing centers. In the presence of collective flow, there is an additional mo- 
mentum transfer orthogonal to the parton trajectory and hence parallel to 
Ht = (— sin ip, cosip). Since the BDMPS transport coefficient denotes the 
squared average momentum transfer per unit pathlength, this contribution 
enters quadratically, |w T (r (£)) ' n T| 2 - For a small collective flow field, this 
quadratic dependence is consistent with the more general ansatz (5.6). 

For an exploratory model study, we use a blast-wave parameterization of the 
hadronic freeze-out stage of the collision [45]. The transverse density distribu- 
tion of the produced matter is specified by 



Here, R x and R y are the extension of the collision region in the reaction 
plane and orthogonal to it, <p s is the azimuthal angle with respect to the re- 
action plane, and we write transverse positions (x, y) in radial coordinates 
x = r cos 0s, y = rsin0 s . We choose a sharp, almost box-like density distribu- 
tion with a s = 0.002. For the flow field u^x), we assume longitudinal Bjorken 
expansion and we use the standard notation 

Ufj,(x) = (cosh rj coshp, sinhp cos 0^, sinhp sin 0b, sinh?? coshp) . (5.12) 

The coordinate 77 denotes the longitudinal space-time rapidity, and we work 
at mid-rapidity i] = 0; 06 defines the orientation which is orthogonal to the 
elliptic freeze-out surface assumed in the model, tan0 s = tan05 (R y /R x ) 2 . 
The transverse flow is parameterized as 




1 



(5.10) 




where f = r(r, <f> s ) denotes a rescaled elliptic position vector, 




(5.11) 



21 



p(r, 4> s ) = f [p + p a cos(20 b )] , 



(5.13) 



where po = 0.88, p a = 0.048 for semi-peripheral Au+Au collisions [45]. For the 
transverse radius parameters, we do not use the extension at freeze-out, but 
the initial transverse radii for an impact parameter b = 7 fm, namely R x = 3.1 
fm and R v = 5.6 fm. 




4 

2 

E 

-2 
-4 



-4 -2 . 2 

x (fm) 



-4 -2 2 4 

x (fm) 









4 












2 

E 






^0 

>< 






-2 






-4 




/ %1 





4 

^2 

E 

>< 

-2 
-4 



x (fm) 



-4 


x (fm) 


4 




L. \ 






1/// i 


5L=4 


-4 


x (fm) 


4 




/ \ 






\ / 




-4 


- 2 , t ° \ 2 
x (fm) 


4 



Fig. 10. Contour plots of the characteristic gluon energy (5.14) as a function 
of the production point ro = (x, y) of the hard parton and for different angles 
<p = 0, vr/4, 7r/2 of its trajectory. The dependence of uj c on the relative flow strength 
If/ Qnf indicates the extent to which hard partons can escape with less energy loss 
on trajectories parallel to the flow field. See text for more details. 
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With this input, we calculate the characteristic gluon energy and average 
transverse momentum squared for a parton trajectory (5.8) in a medium char- 
acterized by its density distribution (5.10) and its collective flow field (5.12). 
With the ansatz (5.9) for the BDMPS transport coefficient, we find 



^ c (ro,V?) 



(qL)(r Q ,<p) = / d£<7(£)0(r (£),£). 
Jo 



(5.14) 
(5.15) 



For a qualitative estimate of the size of parton energy loss, one can use the 
pocket formula AE ~ a s u c [29]. This motivates to investigate u c (r ,(p) as a 
function of the production point r of the hard parton for different orientations 
ip of the parton trajectory. As seen from (5.9), uj c depends linearly on q n f and 
on the relative flow strength qf/q n f- As this flow strength is increased, uo c 
increases for parton trajectories which are not parallel to the flow field. Thus, 
the distortions seen in Fig. 10 provide a first indication of the extent to which 
parton energy loss depends on a transverse flow field and affects the azimuthal 
distribution of inclusive hadron spectra. 

Sernlperipheral AuAu collisions 
v 2 and q„, for R M — 0.5 

0.14 
0.12 
0.1 
0.08 
0.06 
0.04 
0.02 





flow component increases 
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qf/c,, 

Fig. 11. The dependence of elliptic flow V2 and the non-flow component of the 
BDMPS transport coefficient q n f on the relative flow strength qf/qnf, for the case 
of a nuclear modification factor Raa = 0.5 in semi-peripheral Au+Au collisions. 
The calculation is done at fixed transverse momentum px = 7 GeV. 



To estimate the effects of transverse flow, we have calculated from (5.14) and 
(5.15) the relative suppression of hadronic spectra due to medium-induced 
parton energy loss 
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. da med /da vac 

N ^yo^PT) = -^ r / 1 ^r- (5.16) 



K 1 



Our evaluation of (5.16) follows Ref. [9]: we assume a power-law , . 

and we calculate the medium-modification via the quenching weights [9,60,61] 
which depend on (5.14) and (5.15). The integral of (5.16) over r and (p 
weighted with the density of production points determines the nuclear modifi- 
cation factor Raa- We adjust the non-flow component q n f such that Raa = 0.5 
which is the experimentally observed value for semi-peripheral collisions of 
impact parameter b = 7 fm. The results shown in Fig. 11 were obtained for 
Pt = 7 GeV and a s = 1/3. They illustrate two qualitative effects of transverse 
flow: First, low-p^ elliptic flow induces an additional contribution to high-pr 
azimuthal asymmetry. This effect may reduce significantly the discrepancy of 
models of parton energy loss [31,32] in accounting for high-p^ t> 2 . Second, the 
presence of collective flow diminishes strongly the local energy density e oc q^j 
of the medium required for a nuclear modification factor Raa of fixed size. 



6. Conclusion 



In general, hard initially produced partons are not staying in the locally co- 
moving rest frame of the QCD matter generated in a nucleus- nucleus collision. 
Rather, they propagate through a matter which has collective velocity compo- 
nents orthogonal to the parton trajectory. The resulting flow-induced directed 
momentum transfer can modify parton splitting significantly. Here, we have 
studied this effect by calculating the triple-differential medium-induced gluon 
energy distribution (2.9) radiated off a hard parton as a function of gluon 
energy, gluon transverse momentum and azimuthal angle with respect to the 
flow field. Directed momentum transfers lead to a marked asymmetry of the 
medium-induced energy distribution, since partonic fragmentation moves sig- 
nificantly with the direction of the collective flow field, see section 3. 

Based on the medium-induced gluon radiation spectrum (2.3), and simple as- 
sumptions about the dynamical evolution of the matter produced in nucleus- 
nucleus collisions, we have reached several qualitative conclusions of phe- 
nomenological relevance. In particular, as discussed in section 4, flow-induced 
distortions of parton fragmentation will be experimentally accessible in calori- 
metric jet measurements, multiplicity distributions associated to high-pr trig- 
ger particles, and leading two-hadron correlation functions. Moreover, as seen 
in Fig. 11, different combinations of local energy density and collective flow 
can account for the same suppression of single inclusive hadron spectra. This 
illustrates the generic argument of section 5 A that the strength of parton en- 
ergy loss is not governed by the local energy density, but rather by the local 
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energy-momentum tensor (1.2). Flow effects can also contribute appreciably 
to the size of high-p^ ^2 which is underpredicted in recent model compar- 
isons [31,32]. 

In general, the effects of medium-induced parton energy loss depend on time- 
integrated properties of the medium, see section 5 A. Thus, a more quantitative 
study of flow-induced parton energy loss effects requires a realistic model of the 
dynamical evolution of the collision region. It also requires information about 
the spatial distribution of hard processes in the produced matter. Determining 
this information is a challenge which - in an interplay of theory and further 
data analysis at RHIC and the LHC - should come within reach in the near 
future. We hope that our work is of use for further studies in this direction 
and in particular for relating the dynamics of a hydro dynamical medium to 
the dynamics of hard processes in that medium, started in Ref. [62,63]. 
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